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Problem 3 


X 

6 

7 

8 

9 

10 

11 

Y 

0.7782 

0.8451 

0.9031 

0.9542 

1.0000 

1.0414 


Solution : 

b) 

Selective Tow Point Adjacent The Point You Need Fined 


X 0 Xi 


X 

8 

9 

Y 

0.9031 

0.9542 


Fi(x) = f(x 0 ) + f[xi, x 0 ](x-x 0 ) 


f[xi,x 0 ] = 


f(x 0 ) =0.9031 

f(x0)-f(xl) _ 0 . 9542 - 0.9031 
xO-xl 9-8 


f[x 0 ,xi] = 0.0511 


fi(x) = 0.9031+ 0.051(x-8) 
fi(x) = 0.4951 + 0.051 X 
fi(8.3) = 0.4951 + 0.051*(8.3) 
=0.9192 


The Absolute Relative Approximate Error 


0.9191-0.9192 


0.9191 


xlOO 


= 0 . 01 % 


2 





b) 


y. 


X 

7 

8 

9 

Y 

0.8451 

0.9031 

0.9542 


F 2 (x) = f(x 0 ) + f[Xi, x 0 ](x - x 0 ) + f[x 2 , Xi, x 0 ](x - x 0 )(x - Xi). 


f[xi,x 0 ] = 0.0511 

f(x2)-f(xl) f(xl)-f(xO) 0 . 9542 - 0.9031 0 . 9031 - 0.8451 

ffx2 ?( Q1 _ f[x2,xl]-f\xl-x0] _ x 2-xl _ xl-xO _ _ 9-8 8-7 

1 ' ' J x2-x0 x2-x0 9-7 

f[x 2 ,Xi,x 0 ] = -0.0035 

F 2 (8.3) = 0.8451 + 0.0511(x - x 0 ) - 0.0035*(x - x 0 )(x - Xi). 

F 2 (x) = 0.8451 + 0.0511(x - 7) - 0.0035(x - 7)(x - 8) 
f 2 (8.3) = 0.8451 + 0.0511*(8.3-7) - 0.0035(8.3-7)(8.3-8) 

=0.9102 

The Absolute Relative Approximate Error 


=0.97% 


0.9191-0.9102 

0.9191 


x 100 


Problem 7 

The following data defines the sea-level concentration of dissolved oxygen for fresh water as 
function of temperature 


T (c°) 

O (mg/L) 

0 

14.621 

8 

11.843 

16 

9.870 

24 

8.418 

32 

7.305 

40 

6.413 
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Estimate 0(29) using (a) quadratic Lagrange polynomial 

(b) cubic Lagrange polynomial. 


Solution 

a) For quadratic interpolation, the concentration is given by 


0(T) = ^ j L i (T)0(T i ) 


i=0 


= L 0 (T)O( T 0 ) + 4 (T)0(T t ) + L 2 (T)0(T 2 ) 

Since we want to find the concentration at T = 21, and we are using a second order 
polynomial, we need to choose the three data points that are closest to T = 21 that also 
bracket T = 21 to evaluate it. The three points are T 0 =24, 7j = 32, and T 2 = 40. 

Then 

T 0 =24, <9(r 0 ) = 8.418 
7;= 32, 6>(r,) = 7.305 
T 2 = 40, 0(T 2 ) = 6.413 


gives 


T -T. 


j =0 1 0 
j* 0 

' T -I ! N 

' T -T, ' 

T -T 

\ 1 Q 

-u T "‘ 

T -T 

y 1 o 1 2 y 

r 

j 

1=0 1 \ 
i* i 

(t — T 

r j 

f j _ 7y 

T -T 
V J i 

2 T _ 

n 

T -T 

/Vi 2 y 

ii T 

1=0 1 2 
j* 2 

T J 

' t-t \ \ 

7 7 

7 7 


Hence 
0(T) = 


T-T, ^ 

T-T, ' 

O(T 0 ) + 

{ T — T 'l 

' J 0 

T-T, ' 

OlT x ) + 

r t ~ t 0 " 

T-T, ^ 

T -T 

\ 1 Q 1 \J 

T -T 

\ 1 0 1 2 J 

T 

T -T 

1 2j 

T -T 

\ l 2 l 0 J 

T -T 

\ 1 2 


0(T 2 ), T q <T< t 2 
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7X16) 


: (29 - 32X29 - 40) + (29 - 24X29 - 40) 

(24-32)(24-40) (32 - 24)(32 - 40) 

+ (29-24X29-32) 

(40-24)(40-32) 


= (0.2578)(8.418) + (0.8594)(7.305) + (-0.1172)(6.413) 

= 7.694 C 


Problem 7 b : 


T (c°) 

0 (mg/L) 

0 

14.621 

8 

11.843 

16 

9.870 

24 

8.418 

32 

7.305 

40 

6.413 


Solution : 


a) 


0(D = X4(W) 

i =0 

= L 0 C T)O(T 0 ) + L, (7)0(7)) + L 2 (T)0(T 2 ) + L 3 (7)0(7),) 


Then 


gives 


7 0 =16, O(T 0 ) = 9.870 
7; =24, o(7,) = 8.418 
7, =32, o(7 2 ) = 7.305 
7, =40, 0(T 3 ) = 6.413 


3 T — T 

w=n J 


- T -T 

J =0 i 0 1 j 
j* 0 


T -T x ^ 

T -T 1 ' 

' t-t 3 ' 

T -T 

\ 1 Q 1 J 

T -T 

vo 2 y 

T -T 

vo -*3/ 


3 J — J 

h{T) = Y\^jr 

7=0 X ^ j 
M 
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'T 

i 


-T 

- 1 n 


V 


-7, 


oyv 


T 

7) 



T 

T x 




T-Tj 


T -T 

1 2 1 j 


' T 

J2 


-T 

1 n 


V 


-T, 


oyv 


T-T x 

T -T 

. 2 H 


V 


A ^2 




T-Tj 


T -T 

i 3 i j 


{Ti-tJ{t 3 



T 

T 



Hence 


0(T) = 


T-T { 


v 


T n -T 


Vo 

T-T n 


T-T, 


V 


i yv^o 

v 


T n -T 


T-T* 


2 yv-*o 


T n -T. 


O(T 0 ) + 


3 y 


T-T„ 


v 


v r i y v 


r-r. 


v 


T-T, 


T\ ~T 2 ){T{ ~T 3 j 


0(7)) 


\?2 yv 


r-r t 

T -T 
1 2 1 \ 


V 


r-7\ 


yv-^2 7 ) y 


o(t;)+ 


r-r n 


v 


V^3 7; yy 


r-r, 

T -T 

1 3 


v 


r-r. 


yv-^3 7 \j 


0(T,), T 0 <T<T 3 


b) 


0(29) _ (29-24X29-32X29-40) + (29-16)(29-32)(29-40) 


(16 — 24)(16 — 32)(16 — 40) 
(29-16)(29-24)(29-40) 
(32-16)(32-24)(32-40) 
(29-16)(29-24)(29-32) 


(24-16)(24-32)(24-40) 


(7.305) 


(6.314) 

(40-16)(40-24)(40-32) 

= (-0.0537)(9.870) + (0.419)(8.418) + (0.698)(7.305) + (-0.0635)(6.314) 
= 7.695 C 


7.695-7.694 


7.695 


xlOO 


= 0.01299% 
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Problem 11 


X 

-2 

l 

-4 

4 

3 

y 

-l 

2 

-53 

59 

24 


Solution : 


Coordinates 

Ost order 

1-th order 

2-th order 

xo 

f[xo] 



XI 

f[xi] 

f[xo, Xl] 


X2 

f[X2] 

f[xi, X2] 

f[xo, XI, X2] 

X3 

f[X3] 

f[X2, X3] 

f[xi, X2, X3] 


X 

Y 

1 th order 

2 th order 

3 th order 

4 th oi 

-2 

-1 





1 

2 

1 




-4 

-53 

11 

-5 



4 

59 

14 

1 

1 


3 

24 

35 

3 

1 

0 

-1 

4 

5 

6 

1 

0 


From the table is determine the degree of polynomial is 3 


X 

Y 

1 th order 

-2 

-1 


1 

2 

1 


Fi(x) = b G + bi(x-xo) 

F,(0) = -1 +(0 + 2) 

Fi(0) = 1 

From the table at x=0 , y=l 






